Impact of Dresselhaus vs. Rashba spin-orbit coupling on the Holstein polaron 
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We utilize an exact variational numerical procedure to calculate the ground state properties of a 
polaron in the presence of Rashba and linear Dresselhaus spin-orbit coupling. We find that when 
the linear Dresselhaus spin-orbit coupling approaches the Rashba spin-orbit coupling, the Van-Hove 
singularity in the density of states will be shifted away from the bottom of the band and finally 
disappear when the two spin-orbit couplings are tuned to be equal. The effective mass will be 
suppressed; the trend will become more significant for low phonon frequency. The presence of two 
dominant spin-orbit couplings will make it possible to tune the effective mass with more varied 
observables. 



I. INTRODUCTION 

One of the end goals in condensed matter physics is to 
achieve a sufficient understanding of materials fabrica- 
tion and design so as to 'tailor-engineer' specific desired 
properties into a material. Arguably pn-junctions long 
ago represented some of the first steps in this direction; 
nowadays, heterostructures 1 - and mesoscopic geometries 2 - 
represent further progress towards this goal. 

In the field of spintronics, where the spin degree of free- 
dom is specifically exploited for potential applications^ 
spin-orbit coupling^ plays a critical role because con- 
trol of spin will require coupling to the orbital mo- 
tion. Spin orbit coupling, as described by Rashba^ and 
Dresselhaus, 7 is expected to be prominent in two dimen- 
sional systems that lack inversion symmetry, including 
surface states. These different kinds of coupling are in 
principle independently controlled. 8,9 

The coexistence of Rashba and Dresselhaus spin-orbit 
coupling has now been realized in both semiconductor 
quantum wells 4 -" 9 - and more recently in a neutral atomic 
Bose-Einstein condensate. 10 When the Rashba and (lin- 
ear) Dresselhaus spin-orbit coupling strengths are tuned 
to be equal, SU(2) symmetry is predicted to be recov- 
ered and the persistent spin helix state will emergei 4 ' 10 ^ 1 
This symmetry is expected to be robust against spin- 
independent scattering but is broken by the cubic Dres- 
selhaus spin-orbit coupling and other spin-dependent 
scattering which may be tuned to be negligible. 4 

While we focus on the spin-orbit interaction, other in- 
teractions are present. In particular, the electron-phonon 
interaction will be present and may be strong in semi- 
conductor heterostructures. Moreover, optical lattices^ 2 - 
with cold polar molecules may be able to realize a tune- 
able Holstein model J^ The primary purpose of this work 
is to investigate the impact of electron-phonon coupling 
(as modelled by the Holstein modeli 4 -) on the proper- 
ties of the spin-orbit coupled system. We will utilize a 
tight-binding framework; previously it was noted that in 
the presence of Rashba spin-orbit coupling the vicinity 
of a van Hove singularity near the bottom of the elec- 
tron bandit— had a significant impact on the polaronic 
properties of an electron; with additional (linear) Dressel- 



haus spin-orbit coupling the van Hove singularity shifts 
well away from the band bottom, as the two spin-orbit 
couplings acquire equal strength. As we will illustrate 
below, the presence of two separately tunable spin-orbit 
couplings will result in significant controllability of the 
electron effective mass. 



II. MODEL AND METHODOLOGIES 

We use a tight-binding model with dimensionless Hol- 
stein electron-phonon coupling of strength g, and with 
linear Rashba (Vr) and Dresselhaus (Vd) spin-orbit cou- 
pling: 

H = -t J2 ( c l,a c i,» + C ],s c i,s) 

<i,j>,s=ti 

-g^E ^2 clsd^iai+a^+LJE^alai (1) 

where cj s (ci, s ) creates (annihilates) an electron at site 
i with spin index s, and a] (at) creates (annihilates) a 
phonon at site i. The operators Vj, j — 1,2 are written in 
terms of the spin-orbit coupling strengths and the Pauli 
matrices as V\ — Vr& x — Vd<J v , and V<i — Vro v — Vd&x, 
The sum over i is over all sites in the lattice, whereas 
< i,j > signifies that only nearest neighbour hopping 
is included. Other parameters in the problem are the 
phonon frequency, u>e-> an d the hopping parameter t, 
which hereafter is set equal to unity. 

Without the electron-phonon interaction the electronic 
structure is readily obtained by diagonalizing the Hamil- 
tonian in momentum space. With the definitions 



Si = Vr sin(k y ) + V D sin(fc x ), 
S2 = Vr sin(fe x ) + V D sin(fc„), 
we obtain the eigenvalues 

£fc,± = — 2t[cos(k x ) + cos(ky)] ± 2a/S y2 



SI 



(2) 



(3) 



and eigenvectors 



V R /t = 0.5 Vp/t = 0.5 
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FIG. 1. Contour plots for the bare energy bands with Rashba- 
Dresselhaus spin-orbit coupling, for different values of Vr and 
Vd while the sum is kept constant: Vr + Vd = t for these 
cases, (a) Vr = V D = 0.5t, (b) Vr = 0.8i, Vd = 0.2t, (c) 
Vr = 0.9t, Vd = O.li, and (d) Vr = 0.99t, Vd = O.Oli. Note 
the clear progression from a two-fold degenerate ground state 
to a four-fold degenerate one. 
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The ground state energy is 



£ = -*ty/i + (v R + v D y /<&?). 



(4) 



(5) 



Without loss of generality we can consider only Vr > 
and Vd > 0. Either Rashba and Dresselhaus spin-orbit 
coupling independently behave in the same manner, and 
give rise to a four-fold degenerate ground state with 



wave vectors, (k x ,k y ) — (±arctan(-%),± arctan(-^), 
(Vd = 0), and similarly for V D ^ and Vr = 0. With 
both couplings non-zero, however, the degeneracy be- 
comes two-fold, with the ground state wave vectors, 



(k x0 ,k y o) = ±{k ,k ); where k = tan ( 



- U V R + V D 



V2t 



)• 

(6) 

It is clear that the sum of the coupling strengths replaces 
the strength of either in these expressions, so that hence- 
forth in most plots we will vary one of the spin-orbit 
interaction strengths while maintaining their sum to be 
fixed. Similarly, the effective mass, taken along the diag- 
onal, is 



mso 



I 



m VI + (Yr + V D )y(2t 2 ) ■ 



(7) 



where mo = l/(2£) (lattice spacing, a = 1, and K = 1) 
is the bare mass in the absence of spin-orbit interaction, 
and mso is the effective mass due solely to the spin-orbit 
interaction. As detailed in the Appendix, the effective 
mass becomes isotropic when the Rashba and Dressel- 
haus spin-orbit coupling strengths are equal. 

The non-interacting electron density of states (DOS) 
is defined for each band, as 



Ds(e) 



E 

k 



6(e— e ks ) 



(8) 



with s = ±1. 

In FigJ^a) we show the low energy DOS for various 
values of the spin-orbit coupling strengths, Vr and Vd, 
while keeping their sum constant; the low energy van 
Hove singularity disappears for Vr = Vd- Note that only 
D-{e) is shown, as the upper band, with DOS D + {e), 
exists only at higher energies. Furthermore, informa- 
tion concerning the upper band can always be obtained 
through the symmetry 



D + (e) = D_(-e). 



(9) 



In FigJ5][b) we show the value of the density of states 
at the bottom of the band vs. Vd; as derived in the 
Appendix, the DOS value at the minimum energy is given 
by 



D-(E ) = 



1 



1 



2irt 



(Vr+Vd) 2 
2P 



(Vr-Vd) 2 
(Vr + Vd) 2 



(10) 
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FIG. 2. (a)The non-interacting density of states D- (E) 
near the bottom of the band for four values of the spin- 
orbit coupling strengths: (Vr, Vd)/£ = (0.5,0.5) (dot-dashed 
curve), (0.8,0.2) (dotted curve), (0.9,0.1) (dashed curve), 
and (0.99,0.01) (solid curve). Note that for equal coupling 
strengths there is no van Hove singularity at low energies, 
(b) The value of the density of states at the bottom of the 
band (ground state) as a function of Vd (while the total cou- 
pling strength, Vr + Vd, is held constant. The value of the 
density of states achieves a minimum value when Vr = Vd- 
For Vr = or Vd — there is a discontinuity, caused by 
the transition from a doubly degenerate ground state to a 
four-fold degenerate ground state. 



Note that when the coupling strengths are equal, the 
density of states has a minimum. Also note that when 
one kind of spin-orbit coupling vanishes, e.g. Vr = 0, 
or Vd = 0, there will be a discontinuity for the density 
of states (the density of states jumps to twice its value). 
This is caused by a transition from a doubly degenerate 
ground state to a four-fold degenerate ground state. This 
discontinuity will also appear for Vd — or Vr ~ near 
the bottom of the band as can be seen from Figj2^a) for 
Vr = 0.99,V d =0.01. 



III. RESULTS WITH THE 
ELECTRON-PHONON INTERACTION 

As the electron phonon interaction is turned on, the 
ground state energy (effective mass) will decrease (in- 
crease) due to polaron effects. To study the polaron 
problem numerically, we adopt the variational method 
outlined by Trugman and coworkers ; 18 ' 19 which is a con- 
trolled numerical technique to determine polaron prop- 
erties in the thermodynamic limit exactly. This method 
was recently further develope d 20 ' 21 to study the polaron 
problem near the adiabatic limit with Rashba spin-orbit 
coupling^ This case was also studied in Ref. [161 ] using 
the Momentum Average Approximation^ 

In Fig. [31 we show the ground state energy and 
the effective mass correction as a function of the elec- 
tron phonon coupling A = 2<7 2 W£/(47rf)r°- for various 
spin-orbit coupling strengths, but with the sum fixed: 
Vr + Vd = t. These are compared with the results 
from the Rashba-Holstein model with Vd = 0. Here the 
phonon frequency is set to be WE/t = 1.0, which is the 
typical value used in Ref.[16|, and for each value of Vr, 
the ground state energy is compared to the correspond- 
ing result for A = 0. The numerical results are compared 
with results from the MA method and from Lang-Firsov 
strong coupling theor y 23 ' 24 (see Appendix). In Fig. [3Ja), 
the ground state energy crosses over smoothly (at around 
A « 0.8) from the delocalized electron regime to the small 
polaron regime. In the whole regime, the ground state 
energy is shifted up slightly as the Dresselhaus spin-orbit 
coupling, Vd, is increased in lieu of the Rashba spin- 
orbit coupling. We show results for Vd < Vr, as the 
complementary regime is completely symmetric. The 
MA results agree very well with the exact results and 
the Lang-Firsov strong coupling results agree well in the 
A > 1 regime. Similarly, weak coupling perturbation 
theoryil agrees with the exact results for A < 1 (not 
shown). Fig. OJb) shows the effective mass as a func- 
tion of coupling strength; it decreases slightly, for a given 
value of A, by increasing Vd in lieu of Vr. 

All these results are plotted as a function of the elec- 
tron phonon coupling strength, A, as defined above; this 
definition requires the value of the electron density of 
states at the bottom of the band, and we have elected to 
use, for any value of spin-orbit coupling, the value l/(47rf 
appropriate to no spin-orbit coupling. If the actual DOS 
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FIG. 3. (a) Ground state energy difference Egs — Eq vs. A for 
Vn/t — 0.5,0.8, 1.0 and ujE/t = 1.0 while the total coupling 
strength is kept fixed: Vr + Vd — t. Exact numerical results 
are compared with those from the Momentum Average (MA) 
method. Agreement is excellent. Strong coupling results are 
also plotted (in red) by utilizing the Lang-Firsov (LF) strong 
coupling approximation. Agreement in the strong coupling 
regime (A > 1) is excellent, (b) Effective mass m* jmso vs. 
A. MA results are plotted (symbols) with the exact numerical 
results, and again, agreement is excellent. In both (a) and (b) 
the polaronic effects are minimized for Vr — Vd- 



FIG. 4. (a) Ground state energy Egs—Eq as a function of spin 
orbit coupling Vn/t for uiE/t — 0.1, 0.2, 1.0 with weak electron 
phonon coupling, A = 0.32, and moderate spin-orbit coupling, 
Vr + Vd — t. (b) Effective mass m* /mso as a function of 
spin orbit coupling Vb/i for the same parameters. MA results 
are again compared with the exact numerical results, and are 
reasonably accurate for these parameters. 



appropriate to the value of spin-orbit coupling were used 
in the definition of A, then the effective mass, for ex- 
ample, would vary even more with varying Vb vs. Vr 
(see Fig. 2(b)). Moreover, this variation would be more 
pronounced for lower values of uje ■ 

In Fig. 21 we show results for the ground state energy 
and effective mass for different values of the Einstein 
phonon frequency, lue', MA results are also shown for 
comparison. In these plots the electron phonon coupling 
strength is kept fixed and Vd is varied while maintaining 
the total spin-orbit coupling constant. The ground state 
energy has a maximum when the two spin-orbit coupling 
strengths, Vjj and Vr, are tuned to be equal; similarly, 
the effective mass has a minimum when the two are equal. 
As the phonon frequency is reduced the minimum in the 
effective mass becomes more pronounced. The MA re- 
sults track the exact results, and, as found previously, 1 
are slightly less accurate as the phonon frequency be- 
comes much lower than the hopping matrix element, t. 



IV. SUMMARY 

Linear spin-orbit coupling can arise in two varieties; 
taken on their own, they are essentially equivalent, and 
their impact on a single electron, even in the presence of 
electron phonon interactions, will be identical. However, 



with the ability to tune either coupling constant, in both 
solid state and cold atom experiments, one can probe 
the degree of Dresselhaus vs. Rashba spin-orbit coupling 
through the impact on polaronic properties. The primary 
effect of this variation is the electron density of states, 
where the van Hove singularity can be moved as a func- 
tion of chemical potential (i.e. doping) through tuning of 
the spin-orbit parameters. These conclusions are based 
on exact methods (the so-called Trugman method), and 
are not subject to approximations. These results have 
been further corroborated and understood through the 
Momentum Average approximation, and through weak 
and strong coupling perturbation theory. The effect is 
expected to be experimentally relevant since in typi- 
cal materials with large spin-orbit couplings the phonon 
frequency is small when compared to the bandwidth, 
u E /t< 1. 
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Appendix A: Density of States and effective mass 

Expanding Ek,- around the minimum energy Eq, 
by defining k' x = k x ± arctan( Vr ^Z d ), k' y 



k y ± 



arctan( y -"\t y ° ) , we have 

e K - = E + h{K 2 + k' 2 } ± i 2 k' x k' y: 
where 



U =t< 



i , (Vr + Vd) 2 _ (Vr-VdY 

v/l + (Vh + l^) 2 /(2i 2 ) 



(Al) 



(A2) 



and 



U =t< 



(Vr-Vd) 2 
(Vr + V d Y 



y/1 + Wr + VnF/W) 



(A3) 



Note that, with generic spin-orbit coupling, the effective 
mass is in general anisotropic, but when Vd = Vr, it 
becomes isotropic. 

To calculate the density of states at the bottom of the 
band, from the definition, we have 

D_(E + 6E) = ^f dk x f dk y 6(E + 5E-e k ,-), 

where SE is a small amount of energy above the bottom 
of the band, E$. Around the two energy minimum points 
there are two small regions which will contribute to this 
integral. We choose one of them (and then multiply our 
result by a factor of 2) , then use the definitions of k' above 
instead of k, and introduce a small cutoff fc c , which is the 
radius of a small circle around fc m j n . Thus the integral 
becomes 



D_(E +6E) = 2x 



1 



k'dk' 



d6 



1. 



r'l 



8E- {i 1 + -i 2 sm2e}k 



2-Kt / 1 _,_ (Vr + Vd) 2 _ (Vr-Vd) 2 



(A5) 



2/- 



(Vr + V d ) 2 



In the weak electron-phonon coupling regime, perturba- 
tion theory can be applied to evaluate the effective mass; 
the self energy to first order in A is given by 



S woa k(w + iS) = ir\toj E 2J 



1 



k,s=± 



uj + i5 - use - £fe„ 



(A6) 



The effective mass can be obtained through the derivative 
of the self energy 

^^ = 1 - ^-£wcak(^ + iS)\ u =E . (A7) 

m so ouj 

By inserting the expansion of £k,- around the minimum 
energy Eq into Eqn. |A6| and Eqn. [A7J . we obtain the 
effective mass near the adiabatic limit as 



'"-weak 

mso 



X 

1+ 2 



1 



-, , (Vh + Vd) 2 _ {Vr-Vd) 2 
li 2P (Vr + Vo) 2 



(A8) 



The effective mass has a minimum for Vr = Vd while 
Vr + Vd is a constant. 

Appendix B: Strong coupling theory 

To investigate the strong coupling regime of the 
Rashba-Dresselhaus-Holstein model for a single polaron, 
we use the Lang-Firso v 2324 unitary transformation H — 
e s He~ s , where S = g ^2 t n, !CT (aj — a\). Following pro- 
cedures similar to those in Ref. (171), we obtain the first 
order perturbation correction to the energy as 



E 



(i) 

k± 



£fe± 



ff 2 WB, 



(Bl) 



where g is the band narrowing factor, as used in the Hol- 
stein model. To find the second order correction to the 
ground state energy, we proceed as in Ref. (|17l). and find 



E 



(2) 



-4e 



-2s 



2 t 2 + (Vr) 2 + (VdY 
bj E 



x[f(2g 2 )~f(g 2 )]-e- 2 v 2 f(g 2 ) el 



UJ E 



(B2) 



where f(x) = ^ ^ « e*/x[l + 1/x + 2/x 2 + ...]. 

n=l 

Thus the ground state energy, excluding exponentially 
suppressed corrections, is 



E GS = -27rtA(l 



f + (Vr) 2 + (VdY 

(27TtA) 2 



), (B3) 



and there is a correction of order 1/A 2 compared to the 
zcroth order result. Corrections in the dispersion enter 
in strong coupling only with an exponential suppression. 
The ground state energy predicted by strong coupling 
theory has a maximum for Vr — Vd while Vr + Vd is a 
constant. 



